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Abstract 

Theories which have been used to describe the quantized elec- 
tromagnetic field interacting with a nonlinear dielectric medium are 
either phenomenological or derived by quantizing the macroscopic 
Maxwell equations. Here we take a different approach and derive a 
Hamiltonian describing interacting fields from from one which contains 
both field and matter degrees of freedom. The medium is modelled 
as a collection of two-level atoms, and these interact with the elec- 
tromagnetic field. The atoms are grouped into effective spins and the 
the Holstein-Primakoff representation of the spin operators is used to 
expand them in inverse powers of the total spin. When the lowest 
order term of the interaction is combined with the free atomic and 
field Hamiltonians, a Hamiltonian describing a theory of noninter- 
acting polaritons results. When higher order terms are expressed in 
terms of polariton operators standard nonlinear optical interactions 
emerge. These are then compared to the results of phenomenological 
and macroscopic theories. The theory is also used to derive an effec- 
tive Hamiltonian describing counterpropagating modes in a nonlinear 
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medium. 
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1 Introduction 

The study of the propagation of quantized electromagnetic fields in nonlinear 
dielectric media is an area in which there is presently considerable theoretical 
and experimental activity [1,2]. This interest is being driven by both funda- 
mental and practical considerations. On a fundamental level nonlinear optics 
provides an example of an experimentally accessible nonlinear quantum field 
theory, and the light which is produced can exhibit distinctly quantum effects 
such as squeezing and quantum phase diffusion [3-7]. The practical interest 
comes from the proposed use of optical solitons in fibers for long- distance 
communication [8]. This has made necessary a better understanding of the 
role of quantum noise in such systems. 

The first step which must be taken in these problems is the quantization 
of the electromagnetic field in the presence of a nonhnear dielectric. This can 
be done in two ways which we shall call the macroscopic and microscopic ap- 
proaches. In the first, the macroscopic approach, the medium is completely 
described by its linear and nonlinear susceptibilities. No matter degrees of 
freedom appear explicitly in this treatment. One then finds a Lagrangian 
which produces the macroscopic Maxwell equations for the field in a nonlin- 
ear medium. From it one finds the canonical momenta and the Hamiltonian. 
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Quantization is accomplished by imposing the standard equal-time commu- 
tation relations. We explored this approach was in a previous paper [9]. In 
the second approach, the microscopic, a model for the medium is constructed, 
and both the field and matter degrees of freedom appear in the theory. Both 
are quantized. This is what we shall explore here. As we shall see, the result 
is a theory of mixed matter-field modes which are coupled by a nonlinear 
interaction. 

A number of problems arise in the quantized macroscopic theory which 
lead us to examine a microscopic one. The first is the difficulty in including 
dispersion. Dispersion is a result of the fact that the medium is nonlocal 
in time, i. e. the polarization at a given time depends not only on the 
electric field at that time but also on the field at earlier times. It is not 
clear whether nonlocal theories can be encompassed within the canonical 
formalism. (Drummond has shown how it can be done if one is considering 
a narrow frequency interval [10]. He then showed that additional intervals 
could be considered at the expense of introducing additional fields.) Another 
issue is operator ordering. Because the vacuum state of the theory now 
includes the effects of the medium, it is no longer clear that normal ordering 
is the appropriate form to use. Finally, it is expected that the inclusion of 
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losses in the theory will be most easily accomplished at the microscopic level 
[11]. 

The macroscopic and microscopic approaches are complementary; each 
has its own range of application and illuminates the other. The microscopic 
theory we present here, being more fundamental, is more appropriate for 
addressing basic problems such as dispersion and operator ordering. It is, 
however, at the moment restricted to rather simple models of the medium. 
Real media are more easily described by the macroscopic theory in that the 
medium is characterized by only a few parameters. It has been used, for 
example, by Drummond and Carter to study the propagation of quantized 
fields in silica fibers [4] . 

While Drummond and Carter have developed the macroscopic theory 
into a useful tool, much of the work done so far on the nonlinear optics of 
quantized fields has used a phenomenological approach. One assumes that 
the field in a x^^^ medium is described by the nonlinear Schroedinger equation 
and the theory is quantized by imposing equal-space, rather than equal-time, 
commutation relations [5] . A study of the use of equal-space commutations 
has been made by Deutsch [12]. He concluded that for hnear theories they 
lead to the same results as the more standard equal-time relations, but for 
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nonlinear theories it is expected that for some quantities this will no longer 
be the case. These issues are in need of further clarification. 

We shall consider a medium consisting of two-level atoms. This is a model 
which we treated before, but in that case the field was restricted to a single 
mode [13]. This did not allow the inclusion of propagation. We now extend 
that work to include an arbitrary number of modes. This extension is not 
straightforward, because in the single-mode case all of the two-level atoms 
can be treated as one large spin. This is no longer possible when more than 
one mode is present. This necessitates modification of the methods used in 
our first paper. To accomplish this, we divide the medium into little boxes 
whose dimensions are smaller than the shortest wavelength which appears in 
the theory. The atoms in each of these boxes can be treated as a large spin 
and the Holstein- Primakoff representation of the resulting spin operators can 
then be used. This allows us to make a, 1/N expansion of the Hamiltonian 
where N is the number of atoms in a box [14]. A continuum version of the 
theory is obtained by letting the box size go to zero. 

What results is a Hamiltonian containing electromagnetic and matter 
fields. The hnear part of the Hamiltonian is essentially the same as that 
obtained by Hopfield in his treatment of linear dielectric media [15]. It can 
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be diagonalized by transforming from the matter and field modes to mixed 
matter-field modes, known as polaritons. When the nonlinear part of the 
Hamiltonian is expressed in terms of polariton operators the result consists 
of interaction terms which are similar to those which are familiar from non- 
linear optics, i. e. terms describing an intensity-dependent index of refraction 
and third harmonic generation. The effects of dispersion are automatically 
included because the polaritons have a dispersion relation which differs from 
that of photons. The operator ordering in the Hamiltonian is also dictated by 
the theory. The end result is a microscopic theory to which the macroscopic 
and phenonenological theories can be compared. The microscopic theory is 
derived in Sections II-IV and the comparisons with other theories are made 
in Section V. 

2 Hamiltonian for a Medium of Two-Level 
Atoms 

Let us consider a large number, N, of atoms distributed uniformly throughout 
a volume V. These atoms interact with the radiation field whose quantization 
volume is also taken to be V. For simplicity we shall assume that each atom 
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has a single optically active electron. The Hamiltonian describing this system 
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N r 



Here Xj and pj are the position and momentum of the electron on the jth 

— * 

atom which is located at Rj. The potential energy of the electron in the jth 
atom is V{xj — Rj). The charge on the electron is — e where e > 0. Finally, 
we are assuming that we are in the Coulomb gauge so that V • A — 0. 

We now make two approximations. First we wish to specialize to the 
case where only two levels for each atom are included; an upper level \a) 
with energy Eo/2 and a lower level \b) with energy —Eq/2. The projection 



operator onto the two-state space is 



N 



P=\{{\^j){a,\ + h){h,\), 



(2) 



and the effective Hamiltonian for the two-level atoms is obtained by applying 
it to both sides of Eq.(l). Our Hamiltonian now becomes 
"1 



N 

m 



2m 



+ d\[A + (V X Af] } P 



N 



P. 



(3) 



Because we are dealing with optical wavelengths we are also in a position to 
make the dipole approximation and replace Xj, the electron coordinate, by 
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Rj, the position vector of the atom, in the argument of A. This yields 



H = 



N 

m 



1 r --^2 



- I d'r[A + (V X Ay 



N 



p. 



(4) 



Let us now examine the last term of the above equation which represents 
the dipole interaction between the atoms and the field. As a first step we 
expand the vector potential in plane wave modes 



^(a,,e<" + 4_,e-'")^.(£), 



where ur^ = \k\ and the polarization vectors, ex{k), where A 

— * — * 

k ■ ex{k) — and are real. We also adopt the convention 



(5) 

1,2, satisfy 



(6) 



If we now define 



il^xik) = —{a\p\b) ■ ex{k) = ieEo{a\x\b) ■ ex(k), (7) 
m 



then the interaction term becomes 



N 



N 



Y.P[A{R,).p,]P = PJ2J2 



i{^x{k)a^r - ^l{k)a^r))P 



(8) 
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At this point we shall make a number of simplifications. First, we shall 
not write the projection operators explicitly and shall assume that all matter 
operators act only on the two levels, \a) and \b), of each atom. We shall also 
assume that the phases of the atomic wave functions have been chosen so that 
H\{k) is real. Finally, we shall look at the case when only one polarization 
of the field is significantly populated so that the polarization index,A, will be 
dropped. The remaining polarization will be assumed to obey Eq. (6) with 
A = 2. 

We now wish to rephrase the theory in terms of effective spins. To this 
end we divide the total volume, V, into blocks of volume AV. The center 
of the Ith block is located at the position r*;. We suppose that there are Uq 
atoms in each block, where Uq >> 1, and that the dimensions of each block 
are much smaller than an optical wavelength. This means that if we are 
looking at optical phenomena all of the atoms in the Ith block can be treated 
as if they are at r;. The atoms can then be treated as a spin s = no/2 object 
located at fi. Expressing the Hamiltonian in terms of the block variables we 
have 

H-E EoSP + E u;,ay, + ^ E mf + H^m, (9) 
1=1 1=1 
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where 



H^nt = EE + ale-^'-^^)t^^{k){St^ - St\ (10) 

Is) 

We have designated the z component of the spin at fi by Sj; and the raising 
and lowering operators by 5';*^"'"^ and Si~\ respectively. Nj,, which is equal to 
V/AV, is the total number of blocks. 

We can further simplify the A"^ term in the Hamiltonian. If we multiply 
and divide the sum by AV we have 



k 

+ 4ag + a[Vj. (11) 



We have denoted by p the density of atoms which is equal to either N/V or 
no/AV. 

We now have our Hamiltonian in the desired form. The medium is de- 
scribed by spin variables and the field by creation and annihilation operators. 
The next step is to expand it so that we can extract the linear interaction 
and the different orders of nonlinear interaction. We accomplish this with a 
semiclassical, or 1/s, expansion. 
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3 Expansion and Continuum Limit 

— * 

Before proceeding to the actual expansion we must restrict the k summa- 
tions by imposing a high-frequency cutoff. This is necessary so that the 
wavelengths in the theory do not become so small that they violate the con- 
ditions under which the macroscopic theory is valid. In particular, we have 
assumed that the mean spacing between atoms is much less than any wave- 
length in the theory. Therefore, we shall restrict | A; | to be less than ku where 
ku corresponds to a wavelength shorter than those in the optical regime but 
considerably larger than the interatomic separation. 

We would expect that the ground state of the field-atom system would 
be the state in which all of the spins are down, corresponding to all of the 
atoms in their ground states, and no photons present. This is actually the 
case in the semiclassical approximation. A derivation of this fact is presented 
in Appendix A. Here we shall assume this to be the case and shall expand 
our Hamiltonian about the no-photon, all-spin-down state. 

In the case of the spin operators this expansion is implemented by means 
of the Holstein-Primakoff representation [16] in which the spin s operators 
S(^)^ and S^~^ are represented in terms of boson creation and annihila- 
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tion operators, and ( as 



5H = (2.-Ck)n 



(12) 



where [C, C^] — 1- Our convention for the square roots is that when the 
argument is positive so is the square root, and when the argument is negative, 
the square root is i times a positive number. Eqs. (12) then give us the proper 
commutation relations for the spin operators, i. e. 



The excitation number for the boson operators, i. e. the eigenvalue of 
corresponds to S3 + s where S3 is the eigenvalue of S^^\ Therefore, the boson 
vacuum state is the spin state with the spin pointing down. If we are only 
considering states whose excitation number is small then we can expand the 
square roots 



In our model of a nonlinear medium the fraction of atoms in a block which 
is excited will be small because we are off resonance. This corresponds to 
small excitation numbers. Therefore, the use of the expansions in Eq. (14) 



(13) 




(14) 
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is justified. With these approximations our Hamiltonian becomes 

H^Ho + H^^, + H^^„ (15) 

where 

Ho ^ Y: Eo{-s + Cki) + E ^kata^k + i£ 1 d'rA{f)' (16) 
'=1 |fe|<fe,. 

Nb I 



((cVo-c/cf), (18) 

and [0,C/'] = ^w- H^nt represents the linear part of the interaction between 

(2) 

the atoms and the field while H^^l represents the nonlinear interaction. 

We now want to go to a continuum representation of the matter operators. 
Instead of the operators Q and C/ we wish to employ operators, C(r) and C^(r) 
which are functions of a continuous position variable and whose commutation 
relations are 

Note that the operators Q/V AV and (l/VKV have the commutation rela- 
tions 
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^ Sw, (20) 



AV 



which in the hmit — > hmit goes to Eq. (19). This gives us the 
identification 

^0^C(0. (21) 

This relation can be used to find the continuum representation of the 
different terms in the Hamiltonian. For example, the term becomes 

which, after noting that A'^ = 2siVft, gives us 

Ho = ~NEo + EoJ^d'rC\r)m+ E 

\k\<ku 

2 

^'2^^ ^M-£+%-4+°W+4«-.-]- (23) 

\k\<ku 

Similarly, for and i?}^] we find 



= ^ E /^'^K^e*^>4^"*0^A?(C^(r)-C(r1) (24) 

= E / rf'K«.^e*^>4^-^'^-05,KCnr)^C(r) 

-(^^^((r)'), (25) 



where 51^ = jji{k)Jp/2ujk. 
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4 Diagonalization of the Quadratic Part 



So far we have derived a Hamiltonian which contains two interacting fields, 
the electromagnetic field and the matter field. These fields have both linear 
and nonlinear interactions with each other. We have yet to see an interaction 
of the type found in nonlinear optics, that of a field interacting nonlinearly 
with itself. For example, it is not immediately clear how something that 
looks hke the usual description of self-phase modulation would be described 
by our Hamiltonian. 

Something which does look more like the usual theory emerges if we diag- 
onahze the quadratic part of the Hamiltonian, i. e. Hq + H^nt- What results 
is a description in terms of mixed matter-field modes, or polaritons [15]. If 
ij/^j is expressed in terms of polariton operators, it describes a nonlinear in- 
teraction between polaritons. Therefore, instead of a Hamiltonian describing 
the interaction of a photon field with itself, which phenomenological effective 
Hamiltonians in nonlinear optics do, we have a Hamiltonian which describes 
the interaction of a polariton field with itself. 

In order to diagonalize Hq + we introduce the operators 




(26) 
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which imphes that 



C(r) = ^Ee^%- (27) 



These operators have the comutation relations 

[C.-Ci] = %. (28) 

In terms of these operators Hq + h\^^ become 

Ho = -\ne, + e, J2 ClCk+ E 

|fc|<fc„ \k\<ku 

\k\<ku 

and 

= ^ E 9k[HiCl - C-k) + 4^CU - ^k)] ■ (30) 

\k\<ku 

We first note that only the modes k and —k are coupled to each other. 
Therefore, we can diagonalize Hq + h\^^ by considering modes a pair at a 
time and adding up the results. Let be the part of Hq + H\^^ containing 

— * — * 

operators for the modes k and —A;. It can be expressed as 

= P^Biy, (31) 
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where 



V ^-k J 



(32) 



and 



/ 



B 



v 



m 



m 





-m 

-ig^ Eq 



(33) 



In arriving at Eq.(31) we have dropped all constants, such as those arising 
from commutation relations, from the Hamiltonian, and used the fact that 
for our choice of polarization, g_j^ — g^. 

Following Hopfield we want to introduce operators which are linear com- 
binations of the components of u. We shall let a^^ and be these new 
operators where we require 



(34) 



and 



(35) 



19 



Defining the vector 



we have that 



a 



-k 



(36) 



V = Mu, 



(37) 



where M is a 4x4 matrix yet to be determined. The commutation relations 
obeyed by the elements of v and u imply that M has the property that 



M-^ = GM^G, 



(38) 



where the matrix G has the elements Gij — Sij ii i, j — 1,2, and = —Sij 
if i,j — 3,4, with all other matrix elements being zero. Prom Eq. (38) we 
see immediately that M is not unitary. 
Returning to Hj^ we have that 



Hr = v\M-yBM-^v, 



(39) 



which we want to be diagonal. This implies that 



{M-yBM-^ = D 



(40) 
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where D is a. diagonal matrix, or, making use of Eq. (38) and the fact that 

M{GB)M-^ = GD. (41) 

Noting that if D is diagonal so is GD we see that M must diagonalize the 
matrix GB. 

We can reduce finding M to finding the eigenvectors of GB. Let us denote 
the eigenvectors of GB by 77^, where i — 1, . . . , 4, and the unit vector whose 
only nonzero component is the ith by c(i), i. e. c{i)j = Sij. If M satisfies the 
equation 

Mrii = c(z), (42) 

then we find that 

{c{i)\M{GB)M-'\c{j)) = XAj, (43) 

where Aj is the eigenvalue olGB corresponding to r^j. Therefore, if M satisfies 
Eq. (42) it will diagonalize GB. This equation also gives us immediately an 
exphcit representation for 

M-^ = Eh.)(c(^)l- (44) 
It is then possible to use Eq. (38) to find M 

M = G{M-yG = |c(0)(?7,|G (45) 

i=l 
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Finally, let us note that Eq. (38) imposes a condition on the eigenvectors r]i. 
We have from this equation that 

{M-yCM-'^ = G. (46) 

Substitution of Eq.(42) into this result gives us that 

{7^,\Gh) = G,,. (47) 

For i ^ j this condition is automatically satisfied because rji and rjj are 
eigenvalues of GB corresponding to different eigenvalues. For i — j it can be 
imposed as a normalization condition. 

The only task remaining is to find the eigenvectors and eigenvalues of 
GB. Explicit expressions are given in Appendix B along with those for the 
matrix elements of M. Those results along with Eqs. (31) and (37) imply 
that 

H^, = E,{k)a\a^, + ^2(/c)% + E,{k)al^^a_^^ + E2{k)py_^,, (48) 
where 

Ei{k) = l={[El+ujk{uJk + 2Co)\ 

+ [[El - Uk{uk + 2Co)f + 16£;ou;fc^|]V2)V2 (49) 
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- [[E', - UkiiVk + 2Co)]' + leEoUkgl]'/^', (50) 

with Co = {e'^p)/{2mLUk). All of the square roots in the above equations are 
taken to be positive. The energies Ei{k) and E2{k) are just the two branches 
of the polariton energy curve. For large values of ujk we find that Ei{k) — > ujk 
and E2{k) — > Eq. As a;^ — > we find that Ei{k) goes to a value slightly larger 
than 

E,{k)^E, + 2^, (51) 

and £'2 (A;) goes to 

It should be noted that cukCo and uJkQ^ are independent of cuk, and that 
co'fcS'l « Eq and uJkCo « Eq. In Hopfield's analysis the small constant 
term inside the parentheses vanished due to an atomic sum rule. This tells 
us that keeping only one excited atomic level is not a good approximation 
near cuk — 0. Because we are interested in optical phenomena this does not 
present a problem. 

We can finally express Hq + in terms of the polariton operators. 
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Dropping constant terms we have 

^0 + = E {E,{k)a\a^^ + E2{k)plp^). (53) 

\k\<ku 

The effect of the medium appears in two ways in this Hamiltonian. First, 
the operators are mixed matter-field operators, i. e. polariton, not photon, 
operators. Second, the effects of dispersion appear through the energies Ei{k) 
and E2{k) which are not of photon form. Therefore, by including the matter 
degrees of freedom in the theory dispersion emerges naturally, and we avoid 
the problems of trying to quantize a theory which is nonlocal in time. 



5 Nonlinear Interaction 

(2) 

We now turn our attention to the nonlinear part of the Hamiltonian, HlJ^, 
which can be interpreted as describing an interaction between polaritons. 
With the interaction expressed as in Eq. (25), however, this interpretation 
is not obvious. In order to bring it out we begin by expressing it in terms of 
the operators 

^irl = ~^37 E ■ ■ ■ 3k{h+kiM+k2^k^lSl2^ks, 
^ \k\<ku \k3\<ku 

~ h+k2+k3M^kClSk2^k3 + h+k2+k3,ki'^lil2ilsCki 

~ h+k3,kl+k2^lCl^Ck2^kl)- (^^) 
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The next step is to express the photon and atomic operators in terms of 
polariton operators. This leads to a great many terms so it is perhaps best 
to consider a specific physical process and then to select the terms which are 
relevant to it. 

Let us first consider the situation when a single polariton mode, for exam- 
ple, the one corresponding to the operator q;-^^ in the regime where cuk^ > Eq, 

(2) 

is highly excited. The dominant terms in H^^l, at least for times which are 
not too long, will be those in which each of the four operators refers to the 
excited mode. The other terms will have a smaller effect on the time evolu- 
tion because they contain at least one operator for a mode which is initially 
in the vacuum state. Examining the operator transformations in Eq. (37) 
we see that if we only keep the polariton operators and at , then 

Cko ^2i(4)afco C-fco ^ ^Ii(4)4o- (55) 

where have set A — M~^. In the regime we have chosen, > Eq on the 

— * — * — * 

a branch, we find that Aii{ko) is of order one, ^21(^0) and A4i{ko) are of 
order gj^^/ Eq, and A^i^ko) is of order Cq/Eq which is considerably smaller 
than g^^/EQ. Therefore, we shall drop terms containing 7431(^0)- Keeping 
this in mind and making the substitutions indicated in Eq. (55) in Eq. (54) 
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we find 

+ 2\A^^\^Al^a\a-^a\a^^ + {\A^A^Al^ - \A2A''Al^)a^^{alfa-^^ 
+ 2|A4ipA4iKj2(4/], (56) 

— * 

where all of the matrix elements are evaluated at /cq, and we have made use 
of the fact that An^ko) is real and both >l2i(^o) and A4i(/i;o) are imaginary. 

Let us note several things about this expression. First, it looks similar to 
what we would expect an interaction which describes self-phase modulation 
to look like. Therefore, by going to a description in terms of polariton op- 
erators we have recovered a familiar nonlinear optical interaction. Second, 
it is not normally ordered. If we use commutators to bring it into normal 
order we will pick up terms proportional to ctL a^^. These represent small 
shifts to the polariton frequency and can be neglected in most applications. 
Neglecting them gives us 

HlT^ = Xiaifal (57) 



where 



X = -j^An[2{\A2if + 2\A,,\'')A;, + 2{\A,^\^ + 21^21^^41]. (58) 
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The effect of self-phase modulation on a single mode is often described 
by the Hamiltonian 

Hint = \{a)fa\ (59) 

which is superficially similar to Eq. (57). Hoever, in most treatments the 
operators appearing in Eq. (59) are assumed to be photon creation and 
annihilation operators. As is shown by Eq. (57), they should be interpreted 
as polariton operators instead. 

Now let us consider a more complicated situation. Suppose we initially 
have a pulse which is made up of modes on the a branch with wave vectors 

— * 

near k^. In particular, let us assume that all of the wave vectors of the modes 
present in the pulse lie in a small region S about /cq. Keeping only terms in 
h\^^ which contain four excited modes (after transforming to polariton oper- 
ators) we obtain a rather complicated interaction. Each term in it contains 
two creation and two annihilation operators, and four elements of the matrix 

— * 

A{k), each evaluated at a different wave vector. The interaction simphfies 
considerably if we make two approximations. First, we ignore operator or- 
dering, which, as we saw, is tantamount to neglecting small frequency shifts. 
Second, we approximate each matrix element of A by its value at k^. Because 
the spread in wave vectors is small this is a good approximation. With these 
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approximations we find 

^^£r^ = xE E E E h+k.M+k-A'^k'^k.o^ks- (60) 

kes kies k2es k3es 

This situation is often treated by using a phenomenological Hamiltonian 
which describes pulse propagation in a x^^^ medium [5] . The Hamiltonian is 



dx 



(61) 



dx dx 

where the field 4>{x, t) is the field envelope of the pulse in a frame moving at 
the group velocity of the pulse, and it obeys the commutation relations 

[(l)(x,t),(l)''(x',t)]^5(x-x'). (62) 

The equation of motion for the field operator 0(a;, t) resulting from this theory 
is the nonlinear Schroedinger equation. It can also be expressed in terms of 
creation and annihilation operators. Defining the annhilation operator 

a{/3,t)^-^Jdxe'^'(l>{x,t), (63) 

we find for the equation of motion [5] 

.da{p,t) 



i- 



= (5''a{(3, t)+2cj dA J d(32a^{(3i, t)a{(32, t)a{(3 + /3i - (32, t). (64) 



dt 

It should be pointed out that the t in these equations is not really time 
and the x is not really space. The above formalism is derived by an analogy to 
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the classical theory which describes the propagation of a pulse in a nonlinear 
dispersive medium. The pulse is assumed to be propagating in the z direction. 
The variable x is Vgt — where Vg is the group velocity, and t is proportional 
to z. This means that what look like equal-time commutation relations in Eq. 
(62) are actually equal-space commutation relations. The question of when 
a theory using equal-space commutation relations gives the same results as 
a canonically quantized theory has not been fully answered. Deutsch has 
done a preliminary investigation and found that for a linear theory they are 
equivalent [12]. However, for a nonlinear theory he found indications that 
field correlation functions in which the fields are evaluated at different spatial 
points will not be the same in the two theories. Because of the different 
commutation relations, it is rather difficult to directly compare Eq. (64) to 
an equation of motion derived from our interacting polariton theory. 

Instead we shall compare it to one due to Carter and Drummond who 
derived it by quantizing the macroscopic theory [2] . Their theory describes a 
pulse, consisting of modes with wave numbers centered about kc and frequen- 
cies centered about ujc, propagating in a medium of length L in the positive 
z direction. The basic field in this theory is 




(65) 
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which, because of the initial exponential factor, is slowly varying in space 
and time. The annihilation operators are defined in terms, not of the vector 
potnetial and the electric field, but in terms of the dual potential and the 
displacement field [2,13]. They, therefore, imphcitly contain matter degrees 
of freedom and are, in a sense, polariton operators. The field 'if{z,t) obeys 
the equal-time comutation relations 

Mz,t),^^iz',t)]^5iz-z'). (66) 

Like ours, the Carter-Drummond theory describes the pulse in the lab 
frame, i. e. the frame in which the medium is at rest. This removes one 
of the difficulties we had in trying to compare our theory with the phe- 
nomenological nonlinear-Schroedinger equation theory which is formulated 
in a moving frame. Carter and Drummond also reformulate their theory 
in a moving frame in order to make a comparison to the phenomenological 
theory, but we shall use their original lab frame results. In this frame their 
Hamiltonian is 



1 fL 

H ^ - dz 

2 Jo 



dz dz ^ 

' (67) 



where v is the group velocity of the pulse, is proportional to the third 
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order nonlinear susceptibility, and lo" is the second derivative of the frequency 
with respect to the wave number evaluated at kc- The Hamiltonian and the 
commutation relations give, for the equation of motion for the operator afe(i), 

= + v{k - kc) + \uj"{k - kcf]ak 
at I 

- v'^X^ ^ all «fe2 Qik+fei -k2 ■ (68) 

k\ fe2 

It is this equation which we wish to compare to the corresponding equation 
derived from our nonlinear polariton theory. 

Let us find the equation of motion for ot^it). We shall assume that we 

— * 

have a pulse made up of wave vectors near k^ so that the interaction is Eq. 
(60) is appropriate. Because all wave vectors in the pulse are close to AJq, we 
shall expand the polariton energy Ei{k) about /cq 

E,{k) = ^i(^o) + (^^-^o)^ + ^(^^'-(^^-^o)')^ 

+ ^('^^^•^0)^^, (69) 

where 8k = k — kg, ko = ko/kg, and all derivatives are evaluated at kg — \ko\. 
Setting V — dEi/dk, we find 



^ [E^{kg)+vkg-5k + ^^{8p-{5k-kof)v + ]^{5k-ko 



+ 2X E E (70) 



kieSk2es 

31 



If the pulse is one dimensional, i. e. if 5k is always parallel to ko, then Eqs. 
(68) and (70) have the same form. Thus we recover from the microscopic 
theory an equation of motion of the same form as that which arises from the 
quantized macroscopic theory. Note that extending the treatment here to 
broadband pulses is straightforward for microscopic theory; one simply does 

— * 

not expand the polariton energy as a function of k. This extension is more 
complicated for the quantized macroscopic theory. 

We can also examine the more complicated situation in which two modes 
are initially highly excited. Let us look at two cases, one when the modes 
are counterpropagating and one when they are not. The counterpropagating 
case is the more complicated of the two so we shall consider it second. 

Suppose that the two polariton modes on the a branch with wave vectors 

— * — * — * — * 

ka and kb are initially highly excited. Let us assume that ka and kj, are 
both in the x-y plane, and that both modes are polarized in the z direction. 
We shall also assume that ka 7^ —ki,. Keeping only the two excited modes, 
and making the substitutions in Eq. (55) for each mode, the Hamiltonian 
becomes 

Here Xa is simply x with all of the matrix elements evaluated at ka, with a 
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similar definition for Xb, and Xab is given by 

Xab^-^(F-F*), (72) 

where 

F - 4;Wj(l4?r + i4?r)(4f + 4f) + 4f4?(4f + 4i^)] 
4?^4[(l4fr + l4fl^)(4f+4?) 
+4f4f(4f +4?)] (73) 

In the above equation the superscript on the matrix element of A indicates 
whether it is evaluated at ka or ki,. In Eq. (71) one has the usual terms 
which describe cross- and self-phase modulation of the two modes. It should 

— * — * — * 

be noted that if \ka\ — \kb\, then with the above stated conditions on ka and 
kb, and the polarizations, we have that A{ka) = A{kb). This in turn implies 
that Xab — — 4x6, and that the interaction in Eq. (73) can be derived 
from the kind of interaction appearing in the Hamiltonian in Eq. (61) or Eq. 
(67). If \ka\ 7^ \kb\, then this is no longer the case, but will be approximately 

— * — # — * — * 

true if \ka\ — \kb\ is small and the dependence of A(k) on is weak. 

Now let us look at the case of two counterpropagating modes. We shall 

— * — * 

assume that two modes on the a branch with wave vectors kg and —ko, and 
with the same polarization, are initially highly excited. We again want to 
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keep only these two modes in our Hamiltonian, but the situation is now 
considerably more complicated than in our previous cases. This is because 
each of the matter or field operators will now be replaced by a sum of two 
polariton operators, e. g. 

Ck, ^ A2,{ko)aj;^ + A23{ko)al^^, (74) 

where contributions from initially unpopulated modes have been dropped. 
We find that 

+ 4o<«-fc-o] - 2(^21 - ^21^3 + Al,)[{alfal 

where all of the matrix elements are evaluated at kg. In deriving Eq. (75) 
we have made use of the fact that An and ^413 are real, that ^421 and ^23 are 
imaginary (see Appendix B), and that ^13 can be neglected in comparison to 
the other three. We have also dropped terms due to commutators. These are 
similar in form to the terms in Eq. (75) except that they contain only two 
CKjj operators instead of four. They are, therefore, smaller than the terms we 
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have kept by a factor of the order of the number of photons in the and 
q;_^^ modes. If we were to examine the population of initially unpopulated 
modes by the two counterpropagating beams, these terms should be kept. 
In particular, these terms would play a role in the four-wave mixing process 
which is responsible for phase conjugation and squeezing. 

The interaction in Eq. (75) cannot be derived from an effective Hamil- 
tonian like that in Eq. (61). Eq. (75) contains terms with unequal numbers 
of creation and annihilation operators while Eq. (61) does not. This means 
that Eq. (61) cannot be used to treat counterpropagating pulses in a nonlin- 
ear medium, but that a more complicated Hamiltonian, such as that in Eq. 
(75), must be used. 

6 Conclusion 

The justification of the Hamiltonians which are used in the quantum theory 
of nonlinear optics is an important part of placing this theory on a firmer 
basis. We have presented a derivation of an effective Hamiltonian for the in- 
teraction of light and a medium consisting of two-level atoms which contains 
a number of the standard nonlinear optical interactions. This is accomplished 
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by expanding the atomic operators in terms of boson creation and annihila- 
tion operators and diagonahzing the part of the Hamiltonian which describes 
hnear interactions. The result is a theory of interacting polaritons which is 
a nonlinear extension of Hopfield's theory. Because the dispersion relation 
for polaritons is different from that of photons, the effects of dispersion are 
automatically included in this theory. 

The theory presented here provides further justification for the Hamil- 
tonians which emerge from the quantized macroscopic theory. In addition, 
it can be used directly to provide effective Hamiltonians for more compli- 
cated situations than those which have so far been considered. We saw an 
example of this in the case of counterpropagating beams. The microscopic 
theory serves as a useful counterpart to the quantized macroscopic theory in 
providing a description of quantized fields in nonlinear dielectric media. 
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Appendix A 

Here we examine the ground state of the Hamiltonian in the lowest order of 
the semiclassical expansion. The spin operators are replaced by the c-number 
quantities 

^ s cos 61 
^(+) ^ se^'f'^ sin 9i 

Si~^ se-'^'smOi, (76) 

giving us the Hamiltonian 

Hsc = ^os COS 9i+ ^kal^aj^x 

\k\<ku,X 

|fe|<fc„,A '=1 



|fc|<A; 



+ 4a%-. + (-1)Ma«-J' (77) 



where 



2s 
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— * 

We have assumed that fJ,x{k) is real which imphes that D^^^ is also real. In 
principle, we should diagonalize this Hamiltonian which would allow us to 
determine the values of 9i and (pi as well as the field in the ground state 
to lowest order. In practice, we shall be able to obtain what we need after 
completing part of the procedure, and the rest of the diagonalization can be 
completed when additional terms from Hq and i?}^] are included. 

We begin by eliminating the terms linear in the field operators. This can 
be accomplished by shifting the creation and annihilation operators, i. e. by 
setting 

«fcA = hx + ^feA' (79) 

where Zj^^ is a c number, and 

hx = D^kxi-^,:x)HxDkxi-^kxr'- (§0) 

The displacement operator D^^{z) is equal to ex]){zat^ — z*aj^^). We now 
substitute Eq. (79) into Hgc and collect the terms linear in bj^^ and fct^. These 
are 

\k\<k„,X 
2 

+ i^[(-l)'(^fc^A&-fcA + ^-^a^'^a) + 2(4a^£a + Hx^lx) 

+ (-1)'KV-.-a + --.-a4a)]}' (81) 
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where 

L%x-Y.^''''Dc,y (82) 
1=1 

By grouping the terms for k and —k together we find that the expression in 
Eq. (81) will vanish if 

'^'^^^^ + 2m4 + ev) ^ " ^"^^'2m4+% ^^^^ 
Making use of the fact that -^Ij^^ = i~^)^^k\ ^^^^ simplifies to 

Substitution of this result back into Hgc gives us 

Hsc = 51 -^OS COS 9i- Sq^ ^ CUfefct^feg;^ 

+ 4&.-A + (-i)'4^L-J> (85) 



where 



.0= i: (86) 

We can find the semiclassical ground state spin configuration by mini- 
mizing the sum of the first two terms in Eq. (85) with respect to 4>i and 9i. 
That is, we minimize the expression 

J2 EqS COS 9i - Eq. (87) 

1=1 
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We can obtain an estimate of what the minimum configuration will be by 
examining the size of Sq. If £o << NbsEq, then the first term will be the 
dominant one, and the minimum will occur approximately at = tt, for all 
I. This we shall, in fact, find to be the case. 

In order to estimate Eq we first express it in the form 

_ 2{seEof m[{a\x\h) ■ ex{k)f\Ck? (gg) 

where 

Cfe = Ee*''sin^iCos0i. (89) 
1=1 

Taking ku = 2%/ (AVy^^ as the momentum cut off we find that 

E = E sin' Oi cos' 0, < Nl (90) 

This imphes that 

eo< NE^\{a\x\b)\'^m. (91) 

Dividing this by EoNbS, setting the matrix element equal to the Bohr radius, 
and taking £"0 to be the energy of a 500 nm wavelength photon we have that 



^° 5 X 10-3. (92) 



EoNhS 

Therefore, the minimum should be achieved when 9i ^ tt 
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We can, in fact, do better. Let us set 9i — n — 59i and expand the 
expression in Eq. (87) in 59i. We then have that 

EoscosOi ^ -N,Eos + Y.W\ (93) 

1=1 ^ i=\ 

and 

1=1 

<Eos{5x 10-^) Y,W- (94) 
1=1 

From this equation and inequahty it is clear that the minimum occurs at 
66i = 0. Any deviation of 69i from zero causes a larger increase in the first 
term of Eq. (87) than can be compensated for by the decrease in the second. 
This also implies that the shifts Zj^^ are zero and that the operators aj^^ and 
b^^ are the same. The lowest order ground state consists of all the spins 
pointing down and the field in the vacuum state. Givien this information the 
diagonalization of the rest of the Hamiltonian can proceed as in Section 4. 

Appendix B 

— * — * 

The four eigenvalues of the matrix GB are given by Ai = Ei{k), A2 = E2{k), 
A3 = —Ei{k), and A4 = —E2{k). The eigenvector of GB corresponding to 
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the eigenvalue Xj and satisfying Eq. (47) as a normalization condition is 



where 



Vj 



X2j 
\ } 



X4j = 



{El - 



l^^k + Xj\ 

2(a;.|A,|)V2 
{Eo-Xj){u;k + Xjf^' 

LOk — Xj 



nl/2 



{El - X]y + AukEogl 



ujk + Xj 

Eq — Xj 
Eq + Xj 



Xij 



Finally, we can use Eq. (45) to give us that 



(95) 



(96) 



Mij = GiiGjj{rji\c{j)) = GiiGjjX*^, 



(97) 



which, with Eq. (44), imphes that 



(98) 



Let us examine, in particular, the matrix elements which enter into the 
calculation describing two counterpropagating beams (Section 5). We assume 
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that the beams both have wave vectors of magnitude ko and are on the a 



branch. Defining 



we have that 



2\2 



{El - Elf + AuJk,E,gl 



1/2 



(99) 



^11 = 

^13 = 
^23 = 



'^^J^UqEi 

\uko - Ei\d 

"^^^koEi 

£•0 — E\ y E\ 



Eq + EiV ' 



(100) 



and A31 = -Ai3, A33 = -An, A41 = ^23, and A43 = A 



21- 
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